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We present a study of the inclusive photon spectrum from 6.3 million J/ip decays collected 
with the KEDR detector at the VEPP-4M e+e~ collider. We measure the branching 
fraction of the radiative decay J/tp rjc'f, rjc width and mass. Taking into account an 
asymmetric photon line shape we obtain: M,,^ = (2978.1 ± 1.4 ± 2.0) MeV/c^, F^^ = 
(43.5 ± 5.4 ± 15.8) MeV, BiJ/tp rjcj) = (2.59 ± 0.16 ± 0.31)%. 

Keywords: charmonium; radiative decays 

PACS numbers: 13.20.Gd, 13.40.Hq, 14.40.Pq 



1. Introduction 

J/^ — > 77c7 decay is a magnetic dipole radiative transition with photon energy uj 
about 114 MeV, and a relatively large branching fraction of about 2%. This is a 
transition between IS" states of charmonium and its rate can be easily calculated 
in potential models. The transition does not change a spatial part of the wave 
function and its matrix element in the leading approximation equals one. A simple 
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calculation without relativistic corrections gives the resull!-^'i3(J/^/' r]^-f) = 3.05%. 
It was expected that relativistic corrections are of order 20-30%, similarly to the 
case of the electric dipole transitions in charmonium. Therefore, it was surprising 
when in 1986 the Crystal Ball Collaboration measured this branching fraction in 
the inclusive photon spectrum and obtained (1.27 ± 0.36)%P. There are a lot of 
theoretical predictions for this decay rate - - , based on QCD sum rules, lattice 
QCD calculations and so on, but as a rule they give values approximately twice as 
large as the Crystal Ball result. 

This puzzle remained for more than twenty years. Only in 2009 the CLEO 
Collaboration publishecP the result of a new measurement of this branching fraction 
using analysis of 12 exclusive decay modes of 77c. The obtained value B{J/ip 
ijcl) = (1.98 ± 0.09 ± 0.30)% is closer to theory. Combining the Crystal Ball and 
CLEO resuhs, PDQIO got B{ J/^p -J> r]^-/) = (1.7 ± 0.4)% with a scale factor of 1.6. 
In this work we report a result of the new independent measurement. 



2. Photon spectrum 

The photon spectrum in J/ip ^ r]c"f decay can be written aJ^^ 

^ = |a4.^|MPW(.). (1) 

Here M =< ric\joi'-^r/2)\J/'ip > is the matrix element of the transition, jo{x) = 
sin(x)/x, Cc and rric are c-quark charge (in electron charge units) and mass while 
BW{ijf) is a Breit-Wigner function taking into account a nonzero 77c width. The 
matrix element tends to unity when u tends to zero and decreases slowly with the 
photon energy increase. 

CLEO found that the photon line shape of this transition is asymmetric, and a 
Breit-Wigner function alone provides a poor fit to data. Its modification with 
improves the fit around the peak, but gives a diverging tail at higher photon energies. 
To suppress this tail, CLEO used in their fit = exp{--^) with ^ = 65 MeV, 
explaining it by the overlap of the ground state wave functions. However, such a 
form of the matrix element is valid for the wave functions of the harmonic oscillator 
only. In all other potentials, |Afp dependence will be proportional to some negative 
degree of uj when w tends to infinity. 

We tried to fit the CLEO data using another line shape: at photon energy u 
near the resonance the decay probability dV/duj is proportional to uj'^BW{(jj), but 
at higher energies the factor w'^ is replaced with lo. We found that the function 
dT/dbj ^ i^ji^ — yjBW(uj), where loq = — ^^^72 — jg ^jgo suitable. Here the 

' UJUJQ-^-^UJ — CJO) ZJVi , , , ' 

^ ' J/t/' 
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correction factor 7^ ^ is a smooth function near the resonance. Results of 

fits with the CLEO function and our function are shown in Table [TJ One can see 
that results on the ?7c mass, width, and decay rate are close, and the confidence 
level of the fit with our function is also good. Therefore, we use the latter function 
in the analysis of our data. 
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Table 1. Results of the fits to CLEO data using various decay probability functions. 



N^g"^ is the number of si; 


rnal photons in 


the fit. 






dV/duj 


M^^ , MeV/c2 


, MeV 




xVndf (C.L.) 


r^oj^exp{--^)BW(u) 

y2BW{Lu) 


2982.4 ±0.7 
2981.8 ±0.5 


32.5 ±1.8 

33.6 ±1.9 


6142 ± 430 
6494 ± 362 


38.0/38 (0.47) 
39.1/39 (0.47) 



3. KEDR data 

The experiment was performed at the KEDR detectoJ^ of the VEPP-4M 
colhdeM 

The cohider operates with a peak himinosity of about 1.5-10^° cm^^s^^ near the 
J/ip peak energy. Luminosity is measured with single Bremsstrahlung onhne and 
with small-angle Bhabha scattering offline. Two methods of beam energy determi- 
nation are used: resonant depolarization with accuracy of 8 30 keV and IR-light 
Compton backscattering with accuracy ^ 100 kevQ^i, 

This analysis is based on a data sample of 1.52 ± 0.08 pb~^ collected at the J/ip 
peak. Three Jjjt scans were performed. Using a measured beam energy spread as 
well as knowrP^ r( J/?/; — > e+e~) and T{J/ip hadrons), we calculate the J/ip 
production cross section at the peak and get Nj/^ = (6.3 ± 0.3) • 10^. 

Event selection was performed in two steps. At the first step, multihadron decays 
of J/ip were selected with the following cuts: total energy in the calorimeters is 
greater than 0.8 GeV; at least four clusters with the energy greater than 30 MeV 
in the calorimeters are detected; at least one central track in the drift chamber 
is reconstructed; there are no muon tubes activated in the third layer of the muon 
system. These cuts suppress background from the cosmic rays, beam-gas interactions 
and Bhabha events. At the second step, photons in these events were identified. 
A photon is a cluster in the liquid krypton calorimeter which is not associated 
with the reconstructed tracks in the drift chamber and has no time-of-flight (ToF) 
scintillation counters activated in front of it. According to the simulation, the photon 
detection efficiency with the above mentioned cuts is about 28% and is almost 
constant in the investigated range of the photon energies. 



4. Data analysis 

The inclusive photon spectrum and a fit to our data is shown in Fig.[l]. The spectrum 

3 2 

was fit with a sum of the signal having a shape dT/duj ~ ujiuo+{ui-uio)'^ BW{u;), 
convolved with the calorimeter response function (logarithmic normal distributiorP^ 
with aE — 7.4 MeV at 110 MeV and a=-0.33), and background. The background 
shape was taken in the form ln{dN/duj) = a ■ exp{—uj/b) + p2{uj) + c ■ MIP{uj), 
where the first term describes " fake" photons appearing due to nuclear interactions 
of hadrons in the calorimeter and usually having energies less than 60 MeV, the 
second term is the second-order polynomial describing photons arising mainly from 
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Fig. 1. Fit of the inclusive photon spectrum. 

7r° decays, and M IP{uj) is the spectrum of charged particles (a charged particle can 
be also misidentified as a photon). The Breit-Wigner function of the form BW(uj) ^ 
s/((s — M,^J^ + sT^J, where s = Mj^^ — 2ujMj/^, was used in the fit. 

We also tried to fit our data using other line shapes. We used BW{u}) alone, 
uj^BW{uj), the CLEO function, and our function. Again, as first noted by CLEO, 
we found that the BW{uj) function alone gives a shifted value of ijc mass compared 
to other functions: Af(r/c) = 2974.3 ± 1.4 MeV/c^. The function uj^BW{lo) leads to 
a large tail at higher photon energies, giving for the branching fraction B{J/'ijj — >■ 
Vcl) = (7.3 ± 0.5)% (here the decay probability function was integrated up to 
Mj/^/2). The last two functions give close fit results. Since the CLEO function has 
an exponential factor, the result for the branching fraction with this function can 
be considered as a lower limit. The difference between the results obtained with 
the two last functions is used to estimate a systematic error appearing due to the 
unknown line shape. 

5. Systematic errors 

Systematic errors of our measurements are listed in Table [2] 

The uncertainty of the rjc width leads to an error, which we evaluated varying the 
i]c width in the fit by 2.2 MeV (the current PDG error). A systematic error related 
to the background subtraction was estimated by using in the fit a polynomial of the 
third order instead of the second-order one, varying ranges of the fit, and applying or 
not the ToF veto in photon selection. The error in the number of J/ tjj produced was 
evaluated using the known uncertainty of the luminosity measurement. Since the 
cluster multiplicity is different in the simulation and experimental photon spectrum 
in J/ip decays, the error due to the photon detection efficiency was estimated by 
changing by 25% weights of events with small {n < 4) and large {n > 3) track 
multiplicities, and taking different MC generators for the rjc decays. The calibration 
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Table 2. Systematic errors. 



Systematic error 


Mr,^, MeV/c2 


r„e, MeV 


B{J/i>) ^ yVc),% 


Line shape 


0.7 


2.3 


0.15 


rjc width 


0.4 




0.15 


Background subtraction 


0.8 


15.6 


0.17 


Number of J/il> produced 






0.13 


Photon efficiency 






0.08 


Photon energy scale 


1.7 






Total 


2.0 


15.8 


0.31 



of the photon energy scale was performed using 27 decays and tp' — >■ 7Xc j 
'yJ/i' transitions. 

6. Results and conclusions 

A new direct measurement of J/^ — >■ 77c7 decay was performed. We measured 
the r]c mass, width, and branching fraction of J/ip r/c7 decay. The values of 
the branching fraction and rjc mass are sensitive to the line shape of the photon 
spectrum and it should be taken into account during analysis. Our results on the 
77c mass and width are: 

= 2978.1 ± 1.4 ± 2.0 MeV/c^, 
= 43.5 ± 5.4 ± 15.8 MeV. 

Before our experiment these parameters were measured in J/ip and B meson 
decays, as well as in 77 and pp collisions. Mcasurcnncnts of Crystal Ball, MARKS, 
BES, and KEDR were performed using the radiative J/iJj decays, therefore, a mass 
shift should be taken into account. Crystal Ball and KEDR made such a correction, 
but MARKS and BES did not. Therefore we believe that MARKS and BES results 
on the rjc mass should be corrected by approximately 4 MeV towards higher values. 

Our result on the branching fraction of J/tp — rjc"/ decay is 

B{J/ip -> r/c7) = (2.59 ± 0.16 ± 0.31)%. 

It is consistent with that of CLEO, is higher than the old Crystal Ball value and 
close to theoretical predictions. 

The authors are grateful to N. Brambilla and A. Vairo for useful discussions. 
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